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Abstract —This paper studied the fundamental modeling defect 
existing in Poisson-distrihuted cellular networks in which all base 
stations form a homogeneous Poisson point process (PPP) of 
intensity Ab and all users form another independent PPP of 
intensity Ay. The modeling defect, hardly discovered in prior 
works, is the void cell issue that stems from the independence 
between the distributions of users and BSs and “user-centric” cell 
association, and it could give rise to very inaccurate analytical 
results. We showed that the void probability of a cell under 
generalized random cell association is always bounded above 
zero and its theoretical lower bound is exp that can be 

achieved by large association weighting. An accurate expression 
of the void probability of a cell was derived and simulation results 
validated its correctness. We also showed that the associated BSs 
are essentially no longer a PPP such that modeling them as a 
PPP to facilitate the analysis of interference-related performance 
metrics may detach from reality if the BS intensity is not 
significantly large if compared with the user intensity. 

I. Introduction 

Due to the stunning increase of mobile data traffic in cellular 
networks, deploying more and different types of base stations 
(BSs) in cellular networks tends to be an effective means of 
boosting network coverage and capacity. In recent years, mod¬ 
eling and analysis in heterogeneous cellular networks (Het- 
Nets) have gained lots of research attention since traditional 
and simple cellular modeling and analysis cannot correctly 
and accurately characterize the distribution feature of HetNets 
ni-Ei. Stochastic geometry has been shown to be a powerful 
mathematical tool for modeling the random distribution of BSs 
since it makes the mathematical analysis of some network 
performance metrics, such as coverage probability and average 
rate, much more tractable, and it even leads to closed-form 
results in some cases HI 0. 

In the framework of stochastic geometry for cellular net¬ 
works, users are modeled as an independent Poisson point 
process (PPP) in addition to the PPPs formed by all BSs 0 
0. Usually, “user-centric” cell association is adopted in the 
network, i.e., every user tries to associate with its best service 
BS via some cell association schemes. For example, if all users 
are looking for a BS that is able to provide the long-term 
strongest signal power, they will associate with their nearest 
BS if they completely eliminate fading and/or shadowing 
effects on their channels. Such nearest cell association is 
a popular scheme frequently adopted in many prior works 
since the distribution of the distance between a user and its 



Fig. 1. An illustration example of void cells in a cellular network. BSs and 
users are two independent PPPs and the cells of BSs are created by Voronoi 
tessellation. The intensity ratio of users to base stations is 2. 

associated BS can be easily found, and most importantly the 
cells of all macro BSs that consist of the entire network 
plane can be characterized by Voronoi tessellation and thus the 
average number of users in a cell can be accurately estimated 
by either simulated or theoretical computation if the cell of a 
BS is Voronoi-tessellated 0 0. 

A. Motivation and Prior Work 

User-centric cell association has an intrinsic defect that 
gives rise to an unreasonable modeling phenomenon in 
Poisson-distributed cellular networks. Namely, it could make 
BSs not associated by any users in the network, which seldom 
happens in reality especially for macro BSs. This void cell 
(BS) phenomenon is illustrated in Fig. [T] for the case of 
Voronoi-tessellated cells. Users can be imaged to be dropped 
on the plane consisting of the cells so that void cells very likely 
exist in the network if the user intensity is not significantly 
larger than the BS intensity. For example, consider there are 
n users uniformly distributed in a unit-area network which 
is tessellated by n equal-sized cells of BSs, and thus the 
probability of having a void cell is (l — T) that is close 
to e~^ ~ 0.367 as n is very large. This shows that the void 
cell issue cannot be ignored in a cellular network whose BS 
intensity is more or less equal to the user intensity, such as 










hyper-dense small cell networks. 

Despite the important phenomenon of void cells, there are 
very few works that study and model the effect of void cells, 
and almost all prior works on cell association in Poisson 
cellular networks overlook this problem ||3l, |I5], lb), ll9l- lfTT]| . 
Without modeling the effect of void cells, the analyses of 
network performance metrics, e.g. coverage probability and 
average rate, are certainly underestimated since the void BSs 
do not contribute any interference. Although references 112 
m do consider the void cell impact on their models, the 
void probability of a cell they found is only valid for nearest 
cell association, and the associated BSs are still viewed as 
a homogeneous PPP in these two works, which is certainly 
questionable due to the correlation between the associated BSs. 

B. Contributions 

In this paper, our goal is to delve the fundamental corre¬ 
lation between cell association and void cell probability. To 
make the model simple, only a single tier of BSs that form 
a homogeneous PPP of intensity Ab is considered and users 
form another independent PPP of intensity Ay. A random cell 
association (RCA) scheme was proposed to cover several cell 
association schemes, such as nearest cell association, strongest 
received power association, etc.. Our first contribution is to 
theoretically show that the achievable lower bound on the 
void probability of a cell is exp An accurate closed- 

form result of the void probability was derived under the RCA 
scheme for a more practical channel model that characterizes 
path loss as well as the composite effect of Nakagami-m 
fading and log-normal shadowing. According to the derived 
void probability, we can gain some insights on how the void 
probability is affected by channel impairments and how its 
theoretical lower bound is achieved by RCA. Finally, we 
discovered and numerically verified that the associated BSs 
are essentially not a PPP and modeling them as a PPP could 
give rise to fairly inaccurate analysis especially in HetNets and 
small cell networks deployed with high BS intensity. 

II. System Model and Preliminaries 

Consider an infinitely large wireless cellular network on 
the plane in which all base stations form a marked 
homogeneous PPP <I>b of intensity Ab given by 

G {0,1}, iJ, G R+, 

C. cK2,VzGN++|, (1) 

where Bi denotes the iih nearest BS in $b to the origin and 
its location, the cell of Bi is represented by Ct and it contains 
all potential users serviced by BS sfl Vi is a void cell index 
that indicates whether or not at least one user is in the cell of 
Bi, i.e., whether fl <I>u = 0 is true or not - Vi is equal to 
one if Ci n <I>u ^ 0, otherwise zero. Hi is used to characterize 

* The cell of a base station and its statistical properties can be characterized 
once the cell tessellation algorithm (e.g., Voronoi tessellation for a PPP) over 
the plane is designated. 


the downlink channel power gain from Bi to its service user. 
In addition, throughout this paper we assume all cells are 
disjoint, i.e. Ci D Cj =0 for all i ^ j and 
Transmit powers of all BSs are the same and all users form 
another independent homogeneous Poisson point process <I>u 
of intensity Ay. 

To characterize the statistical behavior of downlink channel 
power gains in a much general form, we assume that all HiS 
are i.i.d. random variables whose probability density function 
(pdf) characterizes the composite effect of Nakagami-m fading 
and log-normal shadowing and is given as follows HU 


fnih) = 


T{m)\/2Tra'^ Jo+ 

(In a; — 


-(m-l-l) 


exp 


mh 


2a2 


dec, 


( 2 ) 


where r(ec) = is the gamma function, p and 

are the mean and variance of log-normal shadowing, respec¬ 
tively. Thus, Hi is essentially a gamma-log-normal random 
variable. Note that the mean of Hi is exp (^p + . Without 

loss of generality, the following downlink analysis will be 
based on a typical (reference) user Uq located at the origin. 
Each user associates with a base station in $b by using random 
cell association (RCA). In a given <I>b, the RCA scheme 
associates user Uq with base station Bq by the following rule 


Bo=arg sup {WiH,\\Bi\\ “) , (3) 


where Wi’s are the i.i.d. random weightings for BS association 
and 11 Bill denotes the Euclidean distance between Bi and the 
origin. RCA can cover several BS association schemes. Eor 
example, if the channel state information Hi is available Wi 
can be assigned as Wi = H~^, RCA reduces to the popular 
nearest BS association. If all the weightings of RCA are the 
same constant, RCA is essentially the association scheme of 
strongest received power, which is suitable for non-stationary 
users that cannot acquire the mean received power, e.g., they 
are moving very fast. RCA is a user-centric scheme, i.e., it is 
able to make every user associate with certain base station. In 
other words, no users are blocked out of the network. 

Although user-centric cell association ensures no blocked 
users in the network, it cannot, as we will show later, guarantee 
that every BS is associated with at least one user, i.e., the 
probability that a cell/BS is not associated with any users, 
called the void probability of a cell, is always bounded above 
zero. Since users are a homogeneous PPP of intensity Ay, the 
void probability of a cell in the network can be written as 


P0 = P[E, =0]= [ e-^^"/,(c.)(x)dx, (4) 

where h'(A) is the Lebesgue measure of the bounded Borel 
set Al C and fz{z) is the pdf of random variable Z. In 
order to explicitly calculate p^, one has to know fu{Ci){x) for 
the RCA scheme. The pdf of the cell area of a BS depends 
on the adopted cell association scheme and can be accurately 
characterized by the conservation property of a homogeneous 





PPP introduced in the following subsection. The accurate pdf 
of the cell for RCA will be derived in Section |III] 


A. Random Conservation Property of a PPP 

In this subsection, we introduce the random conservation 
property of a PPP, which specifies how the intensity of a PPP 
is changed after all points of the PPP are transformed by i.i.d. 
random mapping matrices. The random conservation property 
is specified in the following theorem. 


Theorem 1 (Random Conservation Property of a PPP). Sup¬ 
pose $ is an independently marked PPP of intensity measure 
A on , which can be expressed as follows 

$ ^ {{X,,Ti) : X, G G G N++}, (5) 

where Xi denotes node i and its location and Ti : —>■ R'^ 

is the non-singular mapping matrix (operator) of node Xi. For 
all i j, Ti and Tj are (element-wise) i.i.d.. Let be the 
mapped point process on R'^ generated by using the random 
mapping matrix of each point in $, i.e., it is defined as 

$t A {^t A . X, G G G N++}. 

(6) 

Then <1'^ is also a PPP and for any bounded Borel set A C R'^ 
its intensity measure is given by 


At(yf) = A(yf')E 


1 

V'det (TTT) 


(7) 


if ^d'{A') = Ui{A), where A' C R'^ is also a bounded Borel 
set and is the transpose ofT. If^ is homogeneous and has 
an intensity X, is also homogeneous and has the following 
intensity 


AI' = AE 


1 

Vdet(TTT) 


( 8 ) 


Proof: See Appendix lAl 


Remark 1. Theorem\I\is a generalization of the conservation 
property in ( 0 / dm?- In a special case of d = d' = 2, all 
points in are mapped from their corresponding points in a 
homogeneous PPP <i> by scaling them with random variables 
T. In this case, Ti is equal to diag(T, T) and thus A^ is equal 
to AE[ 1 /T 2 ]. 


The random conservation property can significantly reduce 
the complexity of analyzing the statistics of some performance 
metrics induced by a PPP, especially a homogeneous PPP 
with i.i.d. marks. For instance, the RCA scheme in Q can 
be further simplified as 

||i?o*|| ^ ( inf llBlll), (9) 


where = stands for equivalence in distribution, b\ = 
{WiHi)~f Bi, is a homogeneous PPP of intensity Ag = 
AbE {WH)i based on ® in Theorem [T] Hence, for the 
typical user Uq with RCA, its serving distance can be statis¬ 
tically instead found by using the transformed new PPP <i)g. 


In other words, the random conservation property can make 
RCA reduce to nearest BS association, which significantly 
simplifies the analysis of the performance metric of RCA, 
such as the coverage/outage probability of RCA, since many 
existing results of nearest BS association can be applied in this 
context by simply modifying them with an updated intensity 
of the BSs. 


III. Void Probability of a Cell under RCA 


As pointed out in the previous section, any “user-centric” 
association scheme cannot guarantee that there is at least one 
user in each cell, i.e., a void cell could exist in the network. 
This phenomenon can be intuitively interpreted by using a 
Poisson-Dirichlet (Voronoi) tessellation for a PPP of BSs. 
Suppose the Voronoi tessellation is used to determine the cell 
of each BS in $b and all users adopt the nearest BS associ¬ 
ation scheme to connect with their serving BS. This nearest 
associating process can be viewed as the process of dropping 
all users in <i)u on the large plane consisting of the Voronoi- 
tessellated cells formed by $b- Under this circumstance, the 
probability mass function (pmf) of the number of users in a 
cell of $B can be expressed as 


^P[$u(C) =n] =E 


n! 


( 10 ) 


where C denotes a Voronoi cell of a BS in $b, ‘hu(C) 
represents either the user set in cell C or the number of users 
in cell C, and i/(C) is the Lebesgue measure of C. 

Unfortunately, the theoretical result of pn in (doll is un¬ 
known since the pdf of a Voronoi cell area is still an open 
problem Q. However, it can be accurately approximated by 
using a Gamma distribution with some particular parameters 
Q, M- Reference IThl suggests the following Gamma dis¬ 
tribution for 


fv(C){x) 


cX.x 

nox 


(11) 


where r(a;) = is the Gamma function and 

C = I can achieve an accurate pdf of a Voronoi cell area. 
Substituting (fTTT i into (fTOl i yields the following result: 


_ (CAb)^ „+^_i (^Ab+Au)x^ 

- n! r(C) Jo 
i r(n + c) as(cAb)^ 

n\ r(C) (CAb + Au)("+?)’ 


Hence, the pmf of the number of users in a cell for nearest BS 
association has a closed-form expression. The void probability 
of a cell, p 0 , can be found by p„ with the case of n = 0 . 
Namely, 


PH 



(13) 


and this indicates that the intensity of the non-associated BSs 
is AbP 0 that is not negligible especially when Au/Ab is small. 
Most importantly, the results in ([n]l and (Ell are no longer 
accurate for all non-nearest BS association schemes since 


















users do not necessarily associate with their nearest BS, such 
as the RCA scheme mentioned in Section m However, an 
accurate void probability of a BS under the RCA scheme can 
be derived as shown in the following. 

Although the pdf of a Voronori cell is unknown, its mean 
can be shown to be 1/A if the Voronori cell is created by a 
homogeneous PPP of intensity A fT]. As a result, the lower 
bound on in (fTOl i for n = 0 is given by 


P 0 = E 


g-Au--(C) 


> exp 



(14) 


due to Jensen’s inequality. This lower bound on the void 
probability reveals three crucial implications: (i) the void 
probability of a cell is always bounded above zero such that 
there could be a certain number of void BSs in the cellular 
network, (ii) nearest BS association cannot achieve this lower 
bound since its void probability in (foT l is always larger than 
the lower bound, (iii) From an energy-saving point of view, 
the lower bound represents the maximum percentage of all 
consumed energy in the network that can be saved. Later, we 
will theoretically show that, this lower bound can be achieved 
by RCA. 

To derive an accurate void probability of a cell under the 
RCA scheme, we approach this problem from a fundamental 
connectivity point of view and derive the bounds on pg as 
shown in the following theorem. 


Theorem 2. If all users in adopt the RCA scheme defined 
in © to associate with a BS <i>B defined in O, the bounds on 
the void probability of a cell in $b are given by 


1 + 


Au 


where = E \{WH) 


-C 


E 


> P 0 > exp 


Ab 


(15) 


{WH)- 


Proof: See Appendix iBl 


Remark 2. According to the proof of Theorem \2\ the associ¬ 
ation events of different BSs could be correlated such that the 
resulting process of associated BSs is no longer a PPP even 
though its intensity is (1 — P 0 )Ab. 


The lower bound on pg is derived while considering the 
complete independence exists between the non-associated 
events of a BS, whereas the upper bound is obtained by 
approaching the opposite case, i.e., all non-associated events 
of a BS are caused by the farthest user in its cell and thus 
they are highly correlated. Accordingly, it is reasonable to 
conjecture that the upper bound is tightly close to the lower 
bound provided that the cross-correlations between all non- 
associated events are significantly weakened. On the other 
hand, mathematically we have the following 

lim pg = exp (16) 

\ Ab / 

since the upper bound approaches to the lower bound 
exp(—A u/Ab) as (/I goes to infinity and thus the bounds 


in (fTSI l are very tight as becomes large. This intuitively 
reveals that large will “de-correlate” all non-associated 
events, which is an important observation since it indicates 
that the theoretical minimum value of pg is exp and 

it can be achieved by the large ^-fractional moment of WH. 
For example, if W = 1, pg = exp can be achieved 

if channels have high shadowing power. This implies that the 
void probability is reduced when more users join the network 
under a given BS intensity. In other words, when the network 
has a large user population the efficacy of reducing pg by using 
large (jl is apparently undermined such that the performance 
of RCA is similar to that of nearest BS association in this 
case. 

Although the bounds on pg are characterized, an accurate 
result of Pg is still needed since it can help us understand 
how many cells per unit area are void. The following theorem 
renders an accurate closed-form expression of pg under RCA. 

Theorem 3. If RCA is adopted, the void probability of a cell 
can be accurately characterized by 


P0 = 



(17) 


where p = |E \(WH)i E 


Proof: Since (1 -b x/b) ^ > (1 -b x/a) “ for a > b and 
p>C^ we know 


1-b 






-p 


and thus pg in ([TtI i is between the bounds in (fTsT l. Also, for 
nearest BS association the void probability of a BS with C — ^ 
in dnii is very accurate. Therefore, we can conclude that pg 
with p = is an accurate void probability of a BS for the 
RCA scheme since such p reduces to ^ as RCA reduces to 
nearest BS association (i.e., Wi = H~^ for all i). ■ 

According to Remark |2] and Theorem [3 a couple of impor¬ 
tant observations can be summarized as follows. 

> For a cellular network with the condition ^ ^ 1, 
modeling the associated BSs as a homogeneous PPP of 
intensity ( 1 —pg)AB is still acceptable since it just slightly 
looses accuracy in analysis. Hence, a cellular network 
with only macro BSs is hardly affected by the void cell 
problem, but (hyper-dense) small cell and heterogeneous 
networks usually would suffer the void cell problem. 

• Large can considerably reduce the void probability of 
a cell. In other words, if the association weighting Wi is 
designed such that users favor a BS with large channel 
power the void cell problem can be quietly mitigated. 


IV. Simulation Results 

The accurate result of the void probability of a cell has been 
shown in (fTTI l. To verify its correctness, we simulate it for the 
two cases of RCA, i.e., nearest BS association (Wi = 
and strongest received power association (Wi = 1 ) for all i. 














Coverage Probability for Nearest BS Association 



Fig. 2. The void probability of a cell for the two cases of RCA: nearest BS 
association (Wi = H~^) and strongest received power association Wi = 1. 
The path loss exponent a for simulation is 4 and Ay = 370 users/km^. 


For nearest BS association, p in (fTTI i is while for strongest 
received power association we have 


P = 


¥ 



E 



7r(m + ^)r(m 

2 [ r ( TO )]2 



Fig. | 2 ] shows the simulation results of the void probability for 
the two cell association cases. As expected, the void probabil¬ 
ity of the nearest BS association scheme is no longer accurate 
if users associate with their BS with the strongest received 
power. For example, when ^ si 2, the void probabilities 
for the nearest BS association and strongest received power 
association schemes with Rayleigh fading and 8 -dB shadowing 
are 0.2 and 0.14, respectively. The lower bound on the void 
probability is around 0.135. The void probability given in (fTTl i 
indeed accurately coincides with the simulated result and is 
much closer to the lower bound. Also, simulation results and 
(HI both indicate that shadowing can significantly reduce the 
void probability whereas fading just slightly reduces it. 

Next, we simulate how RCA and the void cell problem 
impact the coverage probability of a user. The coverage 
probability of the typical user Uq can be written as 



Fig. 3. The coverage probability of a user for the nearest BS association 
scheme {Wi = l/Hi). The network parameters for simulation are a = 4, 
^ = 0, CF^ = 4dB, /3 = 0.8, m = 1 and Ay = 370 users/km^. 


Coverage Probability for Strongest Received Power Association 



Fig. 4. The coverage probability of a user for the strongest received power 
association scheme {Wi = 1). The network parameters for simulation are 
a = 4, /r = 0, 0-2 = 4dB, /3 = 0.8, m = 1 and Ay = 370 usersdcm^. 


coverage probability is indeed affected by the BS intensity 
since the void probability issue is counted in the model and it 
is a function of the BS intensity. 

V. Concluding Remarks 


Pcou = p [SIRq >/3], (19) 

where SIRy is the signal-to-interference ratio (SIR) evaluated 
at Uq and (3 is the SIR threshold for successful decoding. The 
simulation results for nearest BS and strongest received power 
association schemes are shown in Figs. [3 and @1 respectively. 
As shown in the two figures, the case without considering the 
void cell problem seriously underestimates the real coverage 
probability in the region of ^ < 1. Modeling the associated 
BSs as a PPP of intensity (1—P 0 )Ab can attain a more accurate 
coverage probability, but still not very accurate in the region of 
^ < 1. Note that the coverage probability without considering 
the void cell problem in Fig. [3] does not depend on the BS 
intensity, which has been claimed in Q. However, in fact the 


The fundamental relationship between random cell associa¬ 
tion and void cell probability is exploited in this paper, which 
is an important issue yet overlooked in almost all prior works. 
Random cell association that can cover several cell association 
schemes was proposed, and the accurate result of the void 
probability of a cell is derived and its correctness was verified 
by simulation. Both theoretical and simulation results suggest 
that the void cell problem should be considered and properly 
deal with for cellular networks with small Au/Ab. In order 
to easily characterize the void probability under RCA, only 
one tier of BSs is considered in this paper. Our future work 
will consider an HetNet with more tiers of BSs and focus on 
finding the accurate closed-form result of the void probability 
of a cell in any tier of the HetNet. 

















Appendix 


A. Proof of Theorem |7] 

The following proof essentially generalizes the proof of 
Lemma 1 in insi . First, we define the void probability of ^C\A' 
for any bounded Borel set A' C as the probability that 
is equal to zero where $(,4') denotes the number of 
nodes in $ enclosed in A' . That is. 


= 0 


= exp 



because $ is a PPP. Let Ai = Aii^A') where Ai C is a 
bounded Borel set and the probability that node Xj in is 
enclosed in Ai is 


X] € 


Ai 


= P[T,(X,) e A] = P [T?T,(X,) e Tj{A,)] 
= F[X,g{tJt,)-\A')], 


where (T^Ti)~^(yf') = (T^Ti)~^(dX). By definition, 

is the average of and it can be calculated by 

applying the Campbell theorem as shown in the following; 
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A^jdX^) =E 

Ja 


= E 
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MdX') 


A{dX') 


where (a) follows from the Jacobian determinant of two vol¬ 
umes. Thus, A'l'(yf) = E . ^ A(^') and it follows 


that 


^det(TTT) 


$(^') =0 = exp -E 


v/det(TTT) 
$^(,4) =o' 


-1 


At(dxt) 




since Vd{A) = v^fA'). Since the void probability of a point 
process completely characterizes the statistics of the process, 
is also a PPP 

If $ is homogeneous, then its intensity measure is A(,4') = 
\vdi{A') and thus the result in above can be further simplified 
to A'I'(,4) = AE A/det((T'^T)~i) and it does not 

depend on the location of any node in 4'^. As a result, $Ms 
a homogeneous PPP with the intensity given in ®. 


B. Proof of Theorem ^ 

Suppose each user in cell Ci C R^ of BS Bi has its own cell 
association region which is also a bounded Borel set on R^. 
For example, the cell association region of user Uj is denoted 
by Aj and all cell association regions contain BS Bi. Since 
all users in Ci adopt the RCA scheme, the probability that BS 
Bi is not associated by user Uj G Ci can be expressed as 


P[f.] = 


lim 

u{Aj)^00 


WiHiWUj - Bi\\-^ 


< sup {WiHiWUj - Bi\\ “} , 


where £j denotes the event that user Uj is not associated with 
BS Bi as its association region goes to infinity. According 
to Theorem [T] can be statistically mapped to another 
homogeneous PPP of intensity Ag = ABE[(PFiF)a] 
and 4)g(Alj) is a Poisson random variable with parameter 
Xgiy{Aj). Thus, it follows that 




lim F\\u]-B}\\-°‘< sup 
viAA^oc [ Ble<s>lnAj\B} 



dI > dI* 


where C/J = {WiHi) ^ Uj , is the distance from user UJ 
to BS Bf Uj is the distance from Uj to its nearest BS in 
<i)g. The void probability can be found by calculating the 
following 


P 0 = E < P 


'I’ufCi) 

n 

3=0 


$u(C/) 


Due to the complicate correlation effects between all events 
iSj’s, it is analytically intractable to solve the exact result of 
P 0 . However, the bounds on pg can be characterized in closed- 
form. 

First, the lower bound on p 0 can be found as follows. For 
a given <i)u(C;) we know the following inequality 


‘fu(C0 


n 

$u(C/) 

. 



3>u(Ci) 

> n 


3=0 


which is obtained by considering P[£y|i£i] > F\£j] for 

Ai n Aj f- 0. Since the number of users in Ci is a Poisson 
random variable with parameter i.e., P[$u(C/) = 

n] = g-Aui/(Ci)^ ]-}jg lower bound on pij can be given 

by 


P 0 > E 


SulCi) 

n 


oo 


E 

n =0 



{Xvi3{Ci)r 

gAuu(Ci) 

( 20 ) 


According to the Slivnyak theorem Q, the statistic property 
evaluated at any point in a homogeneous PPP is the same such 




















































that P[iSj] is the same as P[fo] evaluated at user Uq (the origin) 
for all j. Accordingly, (l20l i can be simplified as 


P 0 > exp - P[fo])) 


( 21 ) 


Dl > Dl* 


in which dJ* is the distance 


it* 


and P [£o] = 

from the origin to the nearest BS in $g. The pdf of D, 
fjjuix) = 27r^t^ga;:exp(—7r(^t^ga;;2^ Ijgj^ whereby P[fo] is 
simplified as 

\ 2 '' 


P[£:o] = 1 - E exp Xaioly 


( 22 ) 


Since all users in Ci are uniformly distributed, we can know 
/^t (r) = and thus we can calculate P[fo] in closed- 

form as 


2'KrC^ 

P[£o] = 1 “ y 


dr 


= 1 - 


1 - exp(-C'l'ABi^(Ci)) 


Xbv{Ci) 

Apparently, P[£o] > 1 — Then substituting it into (1^ 

gives us the lower bound in (flST l. 

Next, the upper bound on pg can be obtained as follows. 
The probability that user Uj does not associate with BS Bi is 
upper-bounded by 


P[£j] < lim P 

v{Aj)^oo 


WiHi ini {\\U,-Bi\\-^}< 

Uj 


sup W,H4Uj-B,\\- 

which is acquired by considering the farthest user of BS Bi 
in cell Cl- According to Slivnyak’s theorem and Theorem [T] 
the upper bound on p 0 for a given $(j(C;) can be written as 







p 

n 

$u(Cz) 

< (P 

ot > nt* 

ril ^ 








where rJ denotes the distance from BS Bj to its farthest 
user in C/ = {WH)-^Ci. Let = {U] G C} : 

Uj,Uj G C;} and for given <i>u(Cz) we have 


= p > {Diy 
^ Rh>{Diy 

i=i 



> {Dl 


t*l 2 


^)) 


E[<I.t(ct)] 


'I>u(C,)/Ct 


Au +E[(M^i/)-|]A^g 


(c) 


Ai I 


-CV-tulC,) 


where (a) follows from that fact that = 

1,2,-••} is an one-dimensional Poisson process of inten¬ 
sity AuE[(kLiL)a] and thus is equal to the sum of 

$u(Ci) i-i.d. {WH)~^Rl IfTTlI , (6) is due to first applying 
Jensen’s inequality on the term {WH)~^ and then calculat¬ 
ing the expectation regarding R^, and (c) is obtained from 
(1 — cc)" < (1 — for x G [0,1] and n > 0 and 

E[$u(Cz^)] = ‘l’u(Cz)/C^- Therefore, we have 




n 

$u(Cz) 

_ j=Q 



< 



-ct 


for any given <i)u(Cz) and the upper bound in (fTSl) follows. 
This completes the proof. 
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